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We study linear time fractional diffusion equations in divergence form of time order 
less than one. It is merely assumed that the coefficients are measurable and bounded, and 
that they satisfy a uniform parabolicity condition. As the main result we establish for 
nonnegative weak supersolutions of such problems a weak Harnack inequality with optimal 
critical exponent. The proof relies on new a priori estimates for time fractional problems 
and uses Moser's iteration technique and an abstract lemma of Bombieri and Giusti, the 
latter allowing to avoid the rather technically involved approach via BMO. As applications 
of the weak Harnack inequality we establish the strong maximum principle, continuity of 
weak solutions at t = 0, and a Liouville type theorem. 
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1 Introduction and main result 

Let T > and 51 be a bounded domain in R^. In this paper we are concerned with linear partial 
integro-differential equations of the form 



Here uq — uq{x) is a given initial data for u, A = (aij) is K.^^ ^-valued, Du denotes the spatial 
gradient of u, and stands for the Ricmann-Liouvillc fractional derivation operator with respect 
to time of order a G (0, 1); it is defined by 
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dfv{t,x)—dt / gi^a{t — t)v{t,x) dr, t > 0, x ^ 




where gp denotes the Riemann-Liouville kernel 
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t>0, I3>0. 



*Work partially supported by the European Community's Human Potential Programme [Evolution Equations 
for Deterministic and Stochastic Systems], contract code HPRN-CT-2002-00281, and by the Deutsche Forschungs- 
gemeinschaft (DFG), Bonn, Germany. 



As to applications, equation ((TJ is a special case of problems arising in mathematical physics 
when describing dynamic processes in materials with memory, e.g. in the theory of heat conduc- 
tion with memory, see |24| and the references therein. Time fractional diffusion equations are 
also used to model anomalous diffusion, see e.g. [50]. In this context, equations of the type ([T]) 
are termed subdiffusion equations (the time order a lies in (0, 1); in the case a G (1, 2), which is 
not considered here, one speaks of superdiffusion equations. Time fractional diffusion equations 
of time order a G (0, 1) are closely related to a class of Montroll- Weiss continuous time random 
walk models where the waiting time density behaves as for t -> oo, see e.g. [13], [H], [20] . 

Problems of the type ([T]) are further used to describe diffusion on fractals ([10], [IS]), and they 
also appear in mathematical finance, see e.g. [27] . 

Letting = (0, T) x we will assume that 

(HI) A E LooiflT^R'^''^), and 

N 

\aij{t,x)\'^ < A^, for a.a. {t,x) G fix- 

(H2) There exists > such that 

{A{t, x)^\^) > for a.a. {t, x) € Qt, and all ^ e M^. 

(H3) uo e L2{n). 

We say that a function w is a weak solution (subsolution, supersolution) of (|T]) in Qt, if u belongs 
to the space 

Za:={ve L_^,^{[0,T]; L2in)) n L2i[0,T]; H^{n)) such that 
gi-a*v e C{[0,T];L2{n)), and {gi-a * v)\t=o = 0}, 
and for any nonnegative test function 

r,€H'2\nT) :-ili([0,r];L2(f^))nL2([0,T];iji(17)) (ffi(l]) :^^^^^'(^^)) 
with ri\t=T = there holds 

I lS~ ^t[ffi-a * - Uo)] + {ADu\Di^)'^ dxdt = (<, >) 0. (2) 

Here Lp^ ^ denotes the weak Lp space and /i * /2 the convolution on the positive halfline with 
respect to time, that is (/i * f2){t) = /q fi{t — T)/2(r) dr, t>0. 

Weak solutions of (|T]) in the class Za have been constructed in [36]. Notice also that the 
function uq plays the role of the initial data for u, at least in a weak sense. In case of sufficiently 
smooth functions u and gi_Q * (u — uq) the condition {gi-a * i*)|t=o = implies u\t=o = uq, see 

m- 
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To formulate our main result, let B{x,r) denote the open ball with radius r > centered at 
X £ M^. By UN we mean the Lebesgue measure in M^. For 6 € (0, 1), to > 0, t > 0, and a ball 
B{xo,r), define the boxes 

Q-{to,xo,r) = [to,to + 5Tr'^'") x B{xo,5r), 

Q+ito,xo, r) = {to + (2 - <5)rr2/", io + 2Tr2/") x Bixo,6r). 



Theorem 1.1 Let a € (0, 1), T > 0, and $7 C be a bounded domain. Suppose the assump- 
tions (H1)-(H3) are satisfied. Let further 5 S (0,1), > 1, and t > be fixed. Then for any 
to > and r > with to + 2rr^/" < T, any ball B{xo, rjr) C fi, any < p < 2+jv-,^-2a ' '^'^'^ '^"'2^ 
nonnegative weak supersolution u of in (0, to + 2Tr^/") x B{xo, rjr) with uo > in B{xo, rjr), 
there holds 

I 7 TV / dfiN+i) <C essinf m, (3) 

^llN+l[Q-[to,Xo,r)) jQ_{ta,xo,r) ' g + (to,a;o,r-) 

where the constant C = C(i^, A, (5, r, 77, a, A^, p) . 

Theorem 11.11 states that nonnegative weak supersolutions of ([I} satisfy a weak form of Harnack 
inequality in the sense that we do not have an estimate for the supremum of u on Q^{to, xo, r) 
but only an Lp estimate. We also show that the critical exponent 2+Na-2a optimal, i.e. the 
inequality fails to hold for p > 2+Na-2a ■ 

Theorem 11.11 can be viewed as the time fractional analogue of the corresponding result in 
the classical parabolic case a = 1, see e.g. [121 Theorem 6.18] and [29 . Sending a — 1 in the 
expression for the critical exponent yields 1 + 2/N, which is the well-known critical exponent for 
the heat equation. We would like to point out that the statement of Theorem 11.11 remains valid 
for (appropriately defined) weak supersolutions of ([U on {to, to + 2rr^/") x B{xo, rjr) which are 
nonnegative on (0,io + 2rr^/") x B{xo,rir). Here the global positivity assumption cannot be 
replaced by a local one, as simple examples show, cf. [33]. This significant difference to the case 
a = 1 is due to the non-local nature of d". The same phenomenon is known for integro-differential 
operators like (—A)" with a G (0, 1), see e.g. [TBI . 

As a simple consequence of the weak Harnack inequality we derive the strong maximum 
principle for weak subsolutions of ([Ij, see Theorem 15 . 1 1 below . The weak maximum principle has 
been proven in |32| . even in a more general setting. 

In the classical parabolic case boundedness and the weak (or full) Harnack inequality imply 
an Holder estimate for weak solutions, cf. [5], [TS], [12], [12]. We also refer to [TT] and [H] for 
the elliptic case. In the present situation one cannot argue anymore as in the classical parabolic 
case, due to the global positivity assumption in Theorem ll.il The same problem arises for the 
fractional Laplacian, see |28| . However, in our case it is possible to establish at least continuity at 
t = 0. This is done in Theorem l5.2l in the case uo — 0. It is shown that in this case any bounded 
weak solution u of H]) is continuous at (0,xo) for all xo G and lim(j^')_j(o^^Q) u{t,x) = 0. Thus 
for such weak solutions the initial condition it|t=o = is satisfied in the classical sense. 

As a further consequence of the weak Harnack inequality we obtain a theorem of Liouville 
type, see Corollary 15.11 below. It states that any bounded weak solution of ([T]) on M+ x with 
Uo = vanishes a.e. on M+ x M^. 
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The proof of Theorem 11.11 relies on new a priori estimates for time fractional problems, which 
are derived by means of the fundamental identity (|TOl) (see below) for the regularized fractional 
derivative. It further uses Moser's iteration technique and an elementary but subtle lemma 
of Bombieri and Giusti [5], which allows to avoid the rather technically involved approach via 
BMO-functions. This simplification is already of great significance in the classical parabolic 
case, see Moser [23 and Saloff-Coste [35]. 

One of the technical difficulties in deriving the desired estimates in the weak setting is to 
find an appropriate time regularization of the problem. In the case a = 1 this can be achieved 
by means of Steklov averages in time. In the time fractional case this method does not work 
anymore, since Steklov average operators and convolution do not commute. It turns out that 
instead one can use the Yosida approximation of the fractional derivative, which leads to a 
regularization of the kernel gi—a- This method has already been used in [12] , [30] , [36] , and [32] . 

We point out that the results obtained in this paper can be easily generalized to quasilinear 
equations of the form 

- uo) ~ div a{t,x,u,Du) = b{t,x,u,Du), t G (0,r), x e n, (4) 

with suitable structure conditions on the functions a and b. This is possible, as also known from 
the elliptic and the classical parabolic case, since the test function method used in the proof 
of Theorcni ll.il does not depend so much on the linearity of the differential operator w.r.t. the 
spatial variables but on a certain nonlinear structure, cf. [TT] , [T^] , [IH] ; and [32] . 

In the literature there exist many papers where equations of the type ([T]) , as well as nonlinear 
or abstract variants of them are studied in a strong setting, assuming more smoothness on the 
coefficients and nonlinearities, see e.g. [T], [3], [7], [ID], [H], [21], [33], [31] ■ Concerning the weak 
setting described above one finds only a few results. Existence of weak solutions has been shown 
in [36] in an abstract setting for a more general class of kernels. Boundedness of weak solutions 
has been obtained in [32] in the quasilinear case by means of the De Giorgi technique. With the 
weak Harnack inequality, the present paper establishes a key result towards a De Giorgi-Nash- 
Moser theory for time fractional evolution equations in divergence form of order a S (0, 1). 

We further remark that in the purely time-dependent case, that is for scalar equations of the 
form 

d^{u-uo) +(Tu^O, ie(o,r), 

with cr > 0, a weak Harnack inequality with optimal exponent 1/(1 — a) has been proven in [31j 
for nonnegative supersolutions. Recently, the full Harnack inequality for nonnegative solutions 
has been established in [35]. This, together with the above results, indicates that the full Harnack 
inequality should also hold for nonnegative solutions of ([1]), which is still an open problem, even 
in the case A{t, x) = Id. 

The paper is organized as follows. In Section 2 we collect the basic tools needed for the 
proof of Theorem ll.il These include two abstract lemmas on Moser iterations and the lemma of 
Bombieri and Giusti. We further explain the approximation method for the fractional derivation 
operator and state the fundamental identity (flOl) . which is frequently used in Section 3, where 
we give the proof of the main result. In Section 4 we show that the critical exponent in Theorem 
11.11 is optimal. Finally, Section 5 is devoted to applications of the weak Harnack inequality. 
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2 Preliminaries 



2.1 Moser iterations and an abstract lemma of Bombieri and Giusti 

Throughout this subsection C/o-, < cr < 1, will denote a collection of measurable subsets of 
a fixed finite measure space endowed with a measure /i, such that U^' <Z Ua ii <t' < a. For 
p e (0, oo) and < cr < 1, Lp{Ua) stands for the Lebesgue space Lp{Ua, d^) of all /^-measurable 
functions f:U^^R with |/U^(t/^) := {J^^ \f\P dt,)^/P < oo. 

The following two lemmas are basic to Moser's iteration technique. The arguments in their 
proofs have been repeatedly used in the literature (see e.g. [11], [19], [21], [22], [26], [29]), so it is 
worthwhile to formulate them as lemmas in abstract form, also for future reference. We provide 
proofs for the sake of completeness. 

The first Moser iteration result reads as follows, see also [U Lemma 2.3]. 

Lemma 2.1 Let k > I, p > 1, C > I, and 7 > 0. Suppose f is a fi-measurable function on Ui 
such that 

'C(l 4-/3)'^ 
(cr - a'Y' 

Then there exist constants M — M{C, 7, K,p) and 70 = k) such that 



\f\L,Au^,)<[ ,l^^,' ) \f\L,(u.), 0<a' <a<l, (5>0. (5) 



esssup|/| < ( — — r— - 1 \f\Lp{Ui) for all S G (0,1), p G {0,p]. 



Us 



,(1 -(5)-ro 
Proof: For q> and < cr < 1, let 

*(9,^) = (/ Ifl'dfi)'/'^. 



Let < p < p and S G (0, 1). Set pi — pn^, i = 0,1, . . . and define the sequence {cr^}, i = 0, 1, . . ., 
by o-Q = 1 and 0-^ = 1 — S^-'Xl — S), i = 1,2, . . .; observe that 1 = erg > cri > . . . > cr,; > 

cr,+i > (5 as well as cr^-i — cr^ = 2^'(1 — S), i> 1. Suppose now n S N. By using ([5]) with /? — pi, 
i = 0, 1, . . . , n — 1, we obtain 

HPn,S) <$(p„,CT„) = $(p„_lK,(T„) < ( [2-»(l-^)]7 J" «'(p„-l,fT„-l) 

/ C(2pK"-i)'^ , 
-1 [2-"(1-<5)]J 

^ ^ C(C,p-,7)"«^("'^^ ^--''"'^'^, . . 

=^ i ^3^)^ j HPn-l,^n-l) < ... 

< (c5:?=oO+i)«"^'K7E}ro „ 5)-^5:"r,i«"'y^''$(po,(7o) 
/M(c^y/p 



We let now n tend to 00 and use the fact that 



lim $(p„,(5) = ess sup I/I 
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to get 

esssup|/|<( ) 
Hence the proof is complete. □ 

The second Moser iteration resuh is the fohowing, see also [8, Lemma 2.5]. 

Lemma 2.2 Assume that fJ.(Ui) < 1. Let k > 1, < po < k, and C > 1, 7 > 0. Suppose f is a 
fi-measurable function on Ui such that 

\f\L,uu.,) < ij^^) < a' < C7 < 1, < /3 < ^ < 1. (6) 

Then there exist constants M — M{C, 7, k) and 70 — 70(7, such that 

Proof: Set pi ~ po'* 'i * = 1,2,.... Given 5 G (0,1) we take again the sequence {ct^}, i = 
0, 1,2,.. ., defined by cro = 1 and ai — 1 — 2^-'(l — (5), i > 1. Suppose now n e N. By using 

([6]) with /3 = i = 1, . . . , n, we obtain 

- [2-"(l - ,5)]7«^/P0 [2-("-l)(l - (5)]7'=Vpo *(^'2,cr«-2) < ■■• 
C?^ («+'^' + --- + «") 

Since = pok~'^, we have 

— V = ~ = ^ _ 1) = ^ /J L^ 

Po ~{ Po(k-I) Po(k-I) K-1 Pn PO 



Employing the formula 



^ _l-{n + 1)k" + nK"+i 

j=i ^ ' 



we have further 



+ 1 — = (n + 1) — jK-* 



, - 1 l-(n + 1)k" + 
= (n + 1)k K ^ 



n+l 



K-1 ('t-l)^ 
^n+1 _ + 1)k + n K 

" " (k-I)^ - (k-1)2'' 
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which yields 



Therefore 



1 " K-3 1 1 

l^(n+l-,V<-^(^-l) 



Hpo,s) < 



*(Pn,Cro). 



(1-5)^ 

Given p e (0,po/'^] there exists n >2 such that Pn < P < Pn-i- We then have 

K-1 + -1) 



1 1 k" - 1 



< 



Pn PO PO PO PO 

Pn-1 Po P Pa 



as well as 



<^{pn,cro) ^ $(p„,l) < 
by Holder's inequality and the assumption fi{Ui) < 1. All in all, we obtain 



^{po,6) < 
which proves the lemma. □ 



+ ^) 



HpA), 



The following abstract lemma is due to Bombieri and Giusti [2] ■ For a proof we also refer 
[26l Lemma 2.2.6] and [8] Lemma 2.6] 

Lemma 2.3 Let S, rj G (0, 1), anrf let 7, C 6e positive constants and < (3q < 00. Suppose f 
a positive ^-measurable Junction on Ui which satisfies the following two conditions: 

for all cr, a' , (3 such that Q<5<a'<G<.l and < /3 < min{l, r//3o}- 
(n) 

^x{{\ogf>\})<c^x{Ul)\-^ 



for all A > 0. 
Then 



\f\L,ms)<MKUiy/^'>, 



where M depends only on S, rj, 7, C, and /3q. 
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2.2 The Yoshida approximation of the fractional derivation operator 

Let 0<a<l, l<p< oo, T > 0, and X be a real Banach space. Then the fractional derivation 
operator B defined by 

Bu= ^Cgi_„*7i), D{B)^{u£Lp{[0,T];X): gi^^*ueoH^i{0,T];X)}, 

where the zero means vanishing at t = 0, is known to be m-accretive in Lp( [0, T]; X), cf. [3], 
[6], and [12]. Its Yosida approximations i?„, defined by Bn = nB{n + B)^^, n G N, enjoy the 
property that for any u G D{B), one has BnU Bu in Lp{[0,T]; X) as n — > oo. Further, one 
has the representation 



BnU = - {gi^o.,n*u), u G Lp([0,r];X), n G N, (7) 



d_ 

di 

where gi^a,n = iT-Sa,m ^nd Sa^n IS the unique solution of the scalar-valued Volterra equation 

Sa,nit) + n{Sa,„ * ga){t) = 1, t > 0, n G N, 

see e.g. [30]. Let G Li.;oc(R+) be the resolvent kernel associated with nga, that is 

K,n{i) + n{K^n * ga){t) = nga{t), t>0, neN. (8) 
Convolving ([5]) with gi-a and using g^ * = 1^ we obtain 

] * 9a){t) — n, t > 0, n e N. 

Hence 

gi-a,n^nSa,n^ 9l-a*ha^n, ^ G N. (9) 

The kernels gi-a,n are nonnegative and nonincreasing for all n G N, and they belong to Hl{[0, T]), 
cf. [24] and [30]. Note that for any function / G Lp{[0, T];X), 1 < p < oo, there holds /i^.n*/ ^ / 
in Lp{[0, T];X) as n oo. In fact, setting u = ga * f , we have u G D{B), and 

= ^ (ffl-a.n * ") = 4" (.9l-a * ffa * * /) = h^^n * f Bu = f in ip([0, T]] X) 

at at 
as 71 — ^ oo. In particular, a,n 

in Li([0,r]) as n oo. 
We next state a fundamental identity for integro-differential operators of the form ^(fc * m), 
cf. also _32 . Suppose k G Hl{[Q,T]) and H G C^(R). Then it follows from a straightforward 
computation that for any sufficiently smooth function u on (0,T) one has for a. a. t G (0,T), 

H\u{t))^^ {k * u)it) ^ j^ik* Hiumt) + ( - Hiu{t)) + H'iuit)Ht)) kit) 

{H{u{t - s)) - H{u{t)) - H'{u{t))[u{t - s) - u{t)]^ [-k{s)] ds, (10) 

where k denotes the derivative of k. In particular this identity applies to the Yosida approxi- 
mations of the fractional derivation operator. We remark that an integrated version of (jlOp can 
be found in [TS] Lemma 18.4.1]. Observe that the last term in (|T0|) is nonnegative in case H is 
convex and k is nonincreasing. 

The subsequent two lemmas are also obtained by simple algebra. 
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Lemma 2.4 Let T > and a e (0, 1). Suppose that v e o-ffi^([0,T]) and ip G C^{[0,T]). Then 

{ga*{ipv)){t)^ip{t){ga,*v){t)+ f v{a)d^{gc,{t-a)[ip{t) -ip{a)])da, a.a. te(0,T). 

Jo 

If in addition v is nonnegative and ip is nondecreasing there holds 

(ga * {'fv)){t) > ip{t){ga * v){t) - ga{t - (T)i^{a)v{a)da, a.a.te{0,T). 



Lemma 2.5 LetT>0, ke Hl{[0,T]), v £ Li([0,r]),, and ip e C'^{[0,T]). Then 

^{t) -{k * v){t) = _(fcH<[H)W + k{t^T){^{t)-^{T))v{T)dT, a.a.t£{0,T). 

2.3 An embedding result and a weighted Poincare inequality 

Let r > and be a bounded domain in . For 1 < p < (x we define the space 

Vp ■.= Vpi[0,T]xn)=L2p{[0,T];L2in))nL2i[0,T];H^{n)), (11) 
endowed with the norm 

\u\vp{[0,T]xn) ■= |m|l2p([0,T];L2(O)) + |-Du|L2([o,Tl;L2(n))- 

Set 

.„ 2p + N{p-l) 
^-^P- 2 + N{p~l) 

with = 1 + 2/N. Then Vp ^ L2^{[0,T] x n), and 

l"li2.([o,T]xO) < C{N,p)\u\v^:^[o,T]xn), (13) 

for all u € n L2{[0, T]; i/^ (f^))- This is a consequence of the Gagliardo-Nirenberg and Holder's 
inequality. The case p = oo is contained, e.g., in [18, p. 74 and 75]. The proof given there easily 
extends to the general case. For a more general embedding result (without proof) we also refer 
to [31 Section 2]. 

The following result can be found in [22l Lemma 3], see also [191 Lemma 6.12]. 

Proposition 2.1 Let ip e C(]R^) with non-empty compact support of diameter d and assume 
that < (p < 1. Suppose that the domains {x € M.^ : (p{x) > a} are convex for all a < 1. Then 
for any function u G ^2(1^"'^); 

{u[x) — UipY' (p{x) dx < py-^ — [ \Du{x)\'^ Lp{x) dx , 

R" IV|Li(R") Jr« 

where 

_ Jjjjv u{x)p{x) dx 
/r« dx 
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3 Proof of the main result 



3.1 The regularized weak formulation, time shifts, and scalings 

The following lemma is basic to deriving a priori estimates for weak (sub-/super-) solutions of 
([T]). It provides an equivalent weak formulation of (H]) where the is replaced 

by the more regular kernel gi-a,n (n G N) given in 0. In what follows the kernels /i„ := /iq,„, 
n G N, are defined as in Section [ 



Lemma 3.1 Let a G (0, 1) , T > 0, and fl C be a bounded domain. Suppose the assumptions 
(H1)-(H3) are satisfied. Then u G is a weak solution (subsolution, supers olution) of (QJ) in 
VIt if and only if for any nonnegative function ip G i/jC^) '^'^s ^'^•^ 

(M[ffi-a,« * {u-uo)] + {h„ * [ADu]\DiP)^ dx = (<, >)0, a.a. t G (0,T), n G N. 

For a proof we refer to Lemma 3.1 in 321, where a more general situation is considered with a 
slightly different function space for the solution. The proof of Lemma |3. II is analogous. 

Let u G Za be a weak supersolution of ([T|) in fix and assume that > in f2. Then Lemma 
3.11 and positivity of gi-a,n imply that 

(V9t(5i-a,n *u) + {hn * [ADu] | Z)?/;)) dx > 0, a.a. t G (0,T), n G N, (14) 

for any nonnegative function G H\{^1). 

Let now ti G (0, T) be fixed. For t G {ti,T) we introduce the shifted time s — t — ti and set 
f{s) = f{s + tl), s G (0,r — tl), for functions / defined on {ti,T). From the decomposition 



igi-a,n*u){t,x) ^ / gi^a,nit-T)u{T,x)dT + / gi^a.n{t - t)u{t, x) dr, te(ti,T), 
Jti Jo 

we then deduce that 

dtigi-a.n*u)it,x) ^ ds{gi-a.n*u){s,x) + / gi-a,n{s + ti - t)u{t, x) dr. (15) 

Jo 

Assuming in addition that u > on (0, ti) x fl it follows from p^ . ([T5|) . and the positivity of tp 
and of —gi-a.n that 

(^2pds{gi^a,n*u) + {{hn*[ADu])~\Dij)^dx> 0, a.a. s G (0,T-ti), n G N, (16) 

for any nonnegative function tp G i?2(f^)- This relation will be the starting point for all of the 
estimates below. 

We conclude this section with a remark on the scaling properties of equation Let to,r > 
and xo & K^. Suppose u G Zq, is a weak solution (subsolution, supersolution) of ^ in 
(0,<o^^^") X B{xQ,r). Changing the coordinates according to s = t/r'^^" and y = {x — XQ)/r 
and setting v{s,y) = u{sr'^/°' ,xq + yr), va{y) — uq{xq + yr), and A{s,y) = A{sr'^/°' ,xq + yr), 
the problem for u is transformed to a problem for v in (0, to) x -8(0, 1), namely there holds with 
D = Dy (also in the weak sense) 

dfiv ~ vo) - dW {Ais,y)Dv) ^ i<,>)0, s G (0, to), y e 5(0, 1). (17) 
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3.2 Mean value inequalities 

For cr > we put aB{x, r) :— B{x, ar). Recall that /iat denotes the Lebesgue measure in M^. 

Theorem 3.1 Let a € (0, 1), T > 0, and ft C be a bounded domain. Suppose the assump- 
tions (H1)~(H3) are satisfied. Let r] > and 6 e (0, 1) be fixed. Then for any to e (0,T] and 
r > with to - rjr'^^" > 0, any ball B = B{xo,r) C fi, and any weak supersolution u > e > of 
([IP in (0, to) X B with uq > in B , there holds 

esssup.-i < (^^^±li^)'/V-^k.(r/.), ^ 1, 7 e (0, 1]. 

Here U„ = {to - cn]r^/°',to) x aB, <a <1, C = C{iy,A,S,T],a,N) and tq = To{a,N). 

Proof: We may assume that r ~ 1 and xq = 0. In fact, in the general case we change coordinates 
as i — > t/r^/°' and x ^ (x — xo)/r, thereby transforming the equation to a problem of the same 
type on {0,to/r^/°') x 5(0,1), cf. SectionO 

Fix cr' and a such that S < a' < a < 1 and put Bi ^ aB. For p G (0, 1] we set Vp = Upa- 
Given < p' < p < 1, let ti = <o ^ pcr?? and t2 =to~ p'crr}. Then i) <ti < t2 < to- We introduce 
further the shifted time s = t — ti and set /(s) — /(s + ^i), s £ {O,to — ti), for functions / defined 
on {ti,to). Since uq > in S and m is a positive weak supersolution of ([IJ in (0, to) x B, we have 

(cf. m) 

(vds{gi-a,n *u) + {{hn*[ADu]y\Dv)^ dx > 0, a.a. s G (0,^0 - ii), ?^ e N, (18) 

for any nonnegative function v G Hl{B). For s G {Q,to — ti) we choose the test function 
V = ip^ul^ with /3 < -1 and e C'o(-Bi) so that O<'0<1, ^^ = 1111 p'Bi, suppV' C pBi, and 
< 2/[a{p — p')]. By the fundamental identity (|T0)) applied to fc = 51-a.n and the convex 
function H{y) = -{I + l3)-^y^+P , y > 0, there holds for a.a. {s,x) G (0,to - ti) x B 



i^a,(gi_„,„ * s) > - ^ a.(5i-a,„ * iti+^) + - ^^1+^) 



ds{gi-a,n * u'+n - u'+Pgi-cn. (19) 



l+;3 ^ 1 

We further have 

Dv = 2i;Di;u^ + P^^'^u'^^^ Du. 
Using this and it follows from that for a.a. s G (0,to — ^i) 

1 + P JSi J Bi 

<2 [ {{hni'[ADu]y\ipDiPu'^)dx+ [ 'ijj^u^+^gi-a,,ndx. (20) 

JBi 1 + P 

Next, choose (p G C^([0,to - ti]) such that < < 1, = in [0, (i2 - ^i)/2], (yS = 1 in 
[i2 - ^1,^0 - h], and < < 4/(t2 - ti). Muhiplying 1^ by -(1 + /?) > and by (^(s), and 
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convolving the resulting inequality with yields 

Bi JBi 

<2\1 + I3\ga* [ {{hn*[ADu]y\i^Di^u'')ipdx+\l3\g^* [ i>^u^+P gi^^^^ip dx, (21) 

JBi JBi 

for a. a. s £ (0,to — ^i)- By Lemma 1^?^ 

ga'^{ipds{gi-a,n*bl''^u^'^'^])) dx > / (yj^a * (9s(5i-Q,n * [^/'^M^^'^])) da; 

Bi J Bi 

gc,{s-a)^{a){gi^a,n* [ ^J^u^+^ dx){a) da. (22) 

JBi 

Furthermore, by virtue of 

gi_a,„ * [iP^u'+^] e oHl{[0,to - ti]; Li(Bi)) 

and gi-a,n ~ gi-a * hn as well as ga * gi~a — 1 we have 

* 9.(<7i_o^„ * [V''ui+^]) - * 5i_„,„ * [V'^u'+^D = /in * (^''^1+'^). (23) 

Combining (dH), (|22p . and (1^51) . sending n — cx), and selecting an appropriate subsequence, if 
necessary, we thus obtain 

^ii?u^+^ dx + f3{l + I3)ga* I {ADu\iP^u'^-^Du)ipdx 

Bi Jbi 



< 



2|l + /3|g„*/ {ADu\^D^u^)ipdx+\l3\go,* ijj'^u^+^ gi-a^pdx 



Bi 



+ ga{s- a)tf{a){gi^a* I ij^'^u^^'^ dx){a) da, a.a. s € (0, - ^i)- (24) 

JO J Bi 

Put w = Then Dw = ^^u^^ Du. By assumption (H2), we have 

i3{l + P)ga* [ {ADu\i?u<^-^Du)ipdx >vP{l + I3)ga* ( ^i?u'^-^\Du\^ dx 

JBi JBi 

= T^r9,.* I v^^\Dw\''dx. (25) 

1 + P JBi 

Using (HI) and Young's inequality we may estimate 

2\{ADu\'ipD'ipul^)ip\ < 2Aij\Dip\ \Du\u''^Lp = 2Atp\Dtp\ \Du\u^u^>f 

2 v\l3\ 
- ^P'ip\Dw\'+ ^A^\D^\'ipw\ (26) 



(l + /3)2 ^ ^' ' ly 
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From dM]), (gni), and (HH) we conclude that 



(pii'^w'^ dx + i^'^l^li ga* I (pil/\Dw\^ dx < ga* F, a.a.. s <E {0,to - ti), (27) 

Bi l-L + Pl JBi 



^iPl JSi JBi 



where 



+ 95(s)(gi-Q* / V^"^^ 1^2;) (s) > 0, a.a. s e (0, — ii). 
Jbi 

We may drop the second term in (|27p . which is nonnegative. By Young's inequahty for convohi- 
tions and the properties of (p we then infer that for aU p G (1, 1/(1 — a)) 

(/ mx)w{s,x)]Ux)Pds) < \gc.\LAlo,to~t,]) F{s)ds, (28) 

JBi ' JQ 

where 



l5aU,([o,o-t.]) = r(a)[(a-l)p+l]i/p - r(a)[(a-l)p+l]Vp =^ (29) 

We choose any of these p and fix it. 

Returning to (|27p . we may also drop the first term, convolve the resulting inequality with 
gi-a and evaluate at s = — ti, thereby obtaining 

/*" / i^''\Dw\^dxds< /*" *i^(s)ds. (30) 

Jt2-ti J Bi Jo 

Using 

/ / \D{'iPw)\^dxds<2 I I {^/j^\Dw\^ + \Di}\^w^) dxds 

Jt2~ti J Bi Jt2-ti J Bi 

we infer from (ggi-dSni) that 

+ 4 [ [ \DiP\^w^ dxds. (31) 

Jo Jbi 

We will next estimate the right-hand side of pip . By the assumptions on -0 and ip, and since 
\f3\ > 1, we have 

ptQ — ti /• ^ rto—ti /• 

/ / \D^p\'^w'^ dx ds < / / w^dxds 

Jo Jb, {P-P ) Jo JpBi 
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and 



m < (^^J^^ +|/3|5i-.((t2-ti)/2)) f^^ w'dx 



pBi 

{gi-a* I w'^dx){s), a.a. s G (0, <o — ii)- 



Recall that a > S > 0. So we have 



iD^ dx ds 

pBi 



7 7^ / g2-a{to - ti - t) / w{T,x)'^dxdT 

< G(i^, A, d, 77, aj — ^ — / w dxds. 

[p-pr Jo 



pBi 



Combining these estimates and ([5T|) yields 



IV'w|yp([t2-ti,to-ti]xi3i) < C'l'^, A,(5,77,a,p) "^^^^/^ |w|L2([0,to-ti]xpBi)- 

We apply next the interpolation inequality (jl3p to the lunction ipw and make use of ?/; = 1 in 
p'Bi to deduce that 

kli2.([t2-ti,to-ti]xp'Bi) < C(i/,A,5,?7,a,p, A^) '^^^^/^ |w|L2([o,to-ti]xpBi), (32) 
where the number k > 1 is given in p2p . Since w = and by transforming back to the time 



t, we see that (|32|) is equivalent to 



with (7 = C{v, A, 5, ry, a,p, N). Hence, with 7 = |1 + 

< ) 0<p'<p<l,7>0. 

Employing the first Moser iteration. Lemma [^?T] (with p = 1), it follows that there exist constants 
Mq = MqIv, a, S, rj, a,p, N) and tq = tq(k) such that 

esssupu-i < [ "J \u-^\L^(Vr) for all 9 e (0, 1), 7 e (0, 1]. 



Thus if we take 9 — a' /a and notice that 



1 cr 1 

< 



1-9 a -a' - a -a' 
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we obtain 

esssupu < ( (^_^/)ro J I" \l-,{u.). 7e(0,l]. 
Hence the proof is complete. □ 
We put 

. _ 2 + iVa 

1^ '■— l^l/{l-a) 



2 + Na-2a 



Theorem 3.2 Let a G (0, 1), T > 0, and ft C be a bounded domain. Suppose the assump- 
tions (H1)-(H3) are satisfied. Let rj > and S G (0, 1) be fixed. Then for any to G [0,T) and 
r > with to + rjr"^^" < T, any ball B = B{xQ,r) C fi, any po £ (0,k), and any nonnegative 
weak supersolution u of (Qp in (0,to + rjr^^"') x B with uo > in B, there holds 

^ ^ < < . < 1, < , < po/-.. 

Here = {to, to + arjr^^") x aB, C = C{iy,A,S,7],a,N,po), and To = To{a,N). 

Proof: We proceed similarly as in the previous proof. Without restriction of generality we may 
assume that po > 1 and r = 1. By replacing u with u + e and Uq with Wq + e and eventually 
letting e — 0+ we may further assume that u is bounded away from zero. 

Fix a', a such that 5 < a' < a < 1 and put Bi = aB. For p G (0, 1] we set Vp = U'p„. Given 
Q < p' < p < 1, lei ti — to + p'cTT] and t2 = to + pai], so < to < ti < t2. We shift the time by 
means of s — t ~ to and set f{s) = f{s + to),sG (0, t2 — to), for functions / defined on {to, ^2). 

We then repeat the first steps of the preceding proof, the only difference being that now we 
take /? G (—1,0). Note that, as a consequence of this, simplifies to 

-u'^ds{gi-a,n *u)> - _^ ^ ^ ds{gi-a.n *u^^'^), a.a. {s,x) G (0,t2 - to) X B, 

hence we obtain with -0 G Cg {Bi ) as above 

-TT^ I i^^ds{gi-c,n*u'+^)dx+\l3\ f {{hr. * [ADu])~\^^u^-^Du) dx 

t + P Jbi J Bi 

<2 / {{hn*[ADu])~\i!Di!ul^)dx, a.a. s G (0,^2 -to)- (33) 

Next, choose (/? G C^([0,t2 — ^o]) such that Q < (f < 1, Lp — 1 in [Q,ti — to], (yS = in 
[ti ~to + {t2 - ti)/2, t2 - to], and < -.^ < 4/(^2 - ti). Multiplying dM]) by 1 + /3 > and by 
(p{s), and applying Lemma 12.51 to the first term gives 

ds{gi-a,n*[^'ilJ^u^^'^])dx+]l3]{l+ p) I {ADu]^/j^u^-'^Du)ipdx 

Bi ' JBi 



< 



5i-Q,n(s - cr) (¥'(«)- </'(c^)) ( / i^'^u^^^ dx){a)da 



+ 2(1 + ^) / {ADu]'iPD'iPii'^)(pdx + n„{s), a.a. s G (0,^2 -io), (34) 

Jbi 
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where 

7^„(s) = - |/3|(1 +/3) / {{hr,*[ADu]y - ADu\tl;'^u'^-^Du)ipdx 



Bi 

+ 2(l + /3) / {{hn*[ADu]y - ADu\il}D%l}u^)ipdx, a.a. s e {0,t2 - h). 

JBi 

We set again w — and estimate exactly as in the preceding proof, using (HI), (H3) and 

(Uni), to the resuh 



as(5i-a,n*bV'^w^])rfa;+ / if'>lj'^\Dw\^ dx 

\ J- + P JSi 

< / 9i-a,n{s - a){(p{s) - ip{a)){ 2p^w^ dx) (a) da- 
Jo J Bi 

2A2(l + /3) 



Recall that gi-a,n = gi-a * ^n- Putting 



\D'ip\''Lpw'' dx + nn{s), a.a. s G (0,^2 - <o). (35) 



W{s)^ / ip{s)ip{xfw{s,xydx 

J Bi 

and denoting the right-hand side of pSI) by Fn{s), it follows from ([55]) that 

G„(s) := a,"(/i„ * W){s) + Fn{s) > 0, a.a. s e (0,^2 - to). 
By (j23p and positivity of /i„, we have 

< /l„ * = 5„ * d^ihr, *W)<ga*Gn+ ffa * [-K(s)] + 

a.e. in (0,^2 — ^o), where [y]+ stands for the positive part of y e K. For any p 6 (1, 1/(1 — a)) 
and any G [t2 — tg — [t2 — ti)/4, ^2 — ^o] we thus obtain by Young's inequality 

\hn * W\L^([Q,t.\) < \9a\L^aO,t,]){\Gn\Li([0.t,]) + I h-F„] + |i ^ ([q,*.]) ) . (36) 

Since t* < ^2 ~ < ?7, we have \ga\Lp{[o.t,]) ^ C'i(q;,p, 77) with the same constant as in ([2^. By 
positivity of G„, 



|G„|li([o,u]) = (.9i-a.n * M^)(i*) + / Fnis)ds. 

Jo 

Observe that TZn — >■ in Li([0,i2 — to]) as n — )■ 00. Hence |[— -F'n]+|Li([o.t,]) ^ as n — >• cx). 
Further, 

ngi-a,n(s-cr)((p(s)-(p(cr))( / ip"^ w'^ dx) {a) du ds 
J Bi 

9i^a,n{t* - <j){v{t*) - ^{(y)){ / tjP'uP' dx){(T)da 



'Si 
/,2„,,2 



v{^) / 5i-a,n(s — c) ( / i^P'-uP' dx){a) da ds 

Jo Jbi 

<- I ip{s) j gi-a,n{s - (j){ J ^p'^w'^ dx){a)da ds, 
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since y is nonincreasing. We also know that gi-a,n * — > gi~a * in Li{[0, t2 — io])- Hence we 
can fix some G [^2—^0 — (^2 — ii)/4, ^2—^0] such that for some subsequence {gi-a,nk *W){t^) — >• 
* W){t^:) as fc — > 00. Sending A; — > 00 it follows then from the preceding estimates, 
and from Lp — 1 in [0, ii — ip] that 

"( / [i^{x)w{s,x)Y dxYdsf'" < Ci(a,p,77)((5i-a * W){U) + \F\L,ao^t^-to])) , (37) 
with ^ 

^^^^ ^ 2A_(l_f^ /• |^^|2^y^2^^_^(g)^^^_^ ^ /■ ^2^2 



On the other hand, we can integrate ([55)) over (0,i*) and take the limit as A: — ?■ 00 for the 
same subsequence as before, thereby getting 



^j^\Dw\^dxds< l±l(^{g,_^^W){U) + \F\L,ao,t,-to]))- (38) 



1 + P 

•qj \uw\ ax an ^ 

/O JBi 

Arguing as above (cf. the lines before ([3T|) '). we conclude from ([37)) and ([38)) that 

+ 2(Ci(a,p,77)+ i±^)((5i_„*W^)(<,) + |FU,([o,t,_t„])). (39) 
Since = in [<i - to + (^2 - ^i)/2,t2 - ^o] and t, £ [t2 — to - (<2 — ii)/4,i2 — ^o], we have 



igi-o.*W)iU) <gi.^{{t2~h)/4) I I ifi^w^dxds 

JBi 



r(l-a)(p-p')"(<^'7)" Jo 

Further, 



w 

pBi 



^ dx ds. 



/ / \Dipfw'^dxds< --7T I / w^dxds. 

Jo Jb, (P-P) Jo JpB, 

The term |-F|Li([o.t9-io]) is estimated similarly as in the proof of Theorem 13. II fcf. the lines that 
follow 031])). We obtain 



l + (l + /3) 



w 

pBi 



^ dx ds. 



Notice the additional factor |/3| in the denominator. Combining these estimates we deduce from 
([39)) that 

IV'w|yp([o,ti-to]xSi) < C:{iy,A,S,r],a,p) \w\L2(lo,t2~to]xpBi)- 
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By the interpolation inequality (IT^ and since -0 = 1 in p'Bi, this implies for all /3 G (—1,0) 

kU2.([o,ti-«o]xp'Bi) < C{iy,A,d,rj,a,p,N) -j^j^^^^^ \ML2{[o,t2-to]xpBi), (40) 

where 

2p + jV(p-l) 
2 + iV(p-l) ^ ' ' 

We now fix 1 < p < 1/(1 — a) such that Kp = {pq + k)/2. This is possible because 
1/(1 -a). 

Next, we set 7 = 1 + /3 e (0, 1) and transform back to u to get 



/ C \l/7 

|w|L^jv;,,rfM) < V (p - p/)2 j l"U-7(v^dM): 0<p'<p<l, 0<7< po/k- (41) 



Here, /i = (77Cjjv)^^AiAr+i, wat the volume of the unit ball in R^, and C = C(^, A, (5, 77, a, A'^,po) 
is independent of 7 S (0,po/«^], since |/3| is bounded away from zero. Note that /^(VY) < 1- 

Finally, we employ the second Moser iteration scheme. Lemma 12. 2[ to conclude from (|4T|) 
that there are constants Mq — Mo(iy, A, S, rj, a, N,po) and tq = tq{k) such that 

/ Mo \l/7-l/Po 

|w|L,„(y;,dM) < i (f— ^ j \u\L,iv{,d>^), O<0<1, 0<7< po/«. (42) 

If we take 9 — a' /a and translate back to the measure fJ,N+i, we obtain 

l«U„,.;,, < (M^)'"-"'Vl.,,.;„ 0<7<P./.. (43) 
Since k < k, (|43)) holds in particular for all 7 e (0,po/'*]- This finishes the proof. □ 
3.3 Logarithmic estimates 

Theorem 3.3 Let a e (0, 1), T > 0, and ft C be a bounded domain. Suppose the assump- 
tions (H1)-(H3) are satisfied. Let r > and 6, rj £ (0, 1) be fixed. Then for any tg > and 
r > with to + rr^/" < T, any ball B = B(xo,r) C fl, and any weak supersolution u > e > of 
OP m (0,io + -rr2/°) x B with uq > in B , there is a constant c = c{u) such that 

HN+i{{{t,x)eK^:logu{t,x)>c + X})<Cr^^°'fiN{B)X-\ A > 0, (44) 

and 

fiN+i{{it,x) e K+ ■.logu{t,x) < c~ X}) <Cr^^°'^lNiB)X-\ A > 0, (45) 

where K- = (to, to + r]Tr'^/°') x SB and — (to + r^rr^/", to + rr^^") x SB. Here the constant 
C depends only on S,ri,T, N,a,i>, and A. 

Proof: Since uq > in B and u is a positive weak supersolution we may assume without loss of 
generality that uq = and to = 0. In fact, in the case to > we shift the time as t — > t — to, 
thereby obtaining an inequality of the same type on the time-interval J :— [0,rr^/"]. Observe 
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that the property gi^a * u G C([0,<o + Tr^^"]; L2{B)) miphcs gi-a * u & C{J; L2{B)) for the 
shifted function u(s, x) = u{s + to, x). So we have 

J (ydt{gi^a,n *u) + (hn * [ADu]\Dv)^ dx > 0, a.a. teJ,neN, (46) 

for any nonnegative test function v 6 H2{B). 

For t € J we choose the test function v — ip'^u"^ with -0 G Cl{B) such that suppt/' C -B, 
-i/j = 1 in 5B, < V' ^ Ij l-DV'l 2/[(l — 5)r\ and the domains {x ^ B : ip{x)'^ > b} are convex 
for all 6 < 1. We have 

Dv 2ipDipu~'^ - ij/u-^Du, 
so that by substitution into we obtain for a.a. t <E J 

'ij?u^^dt{gi-a,n*u)dx+ i {ADu\u^^Du)^p'^ dx 

<2 {ADu\u-^ ip Dip) dx + Unit), (47) 
Jb 

where 

7e„(t) = / [h„* [ADu] - ADu\Dv) dx. 
Jb 

By (HI) and Young's inequality, 

\2{ADu\u-^ijDi;) \ < 2A'iIj\Di/j\ \Du\u-^ < ^ iP^\Du\^u-^ + ^ A^jDi/'p. 
Using this, (H2) and \DiIj\ < 2/[(l - 6)r], we infer from (gTl) that for a.a. t e J 

- I ^'u~%{g,_^^^*u)dx+'^ I \Du\^u-^^^dx< ^.f^^lfl +7^„(0. (48) 

Setting w — log it we have Dw — u^^Du. The weighted Poincare inequality of Proposition 12.11 
with weight ip'^ yields 

I {w- Wfi^^dx < f \Dw\''i;^dx, a.a. t G J, (49) 

Jb . r W dx Jb 



where 



. N fn wU, x)ih(x)'^dx 
^(^) = r //l2\; ' a.a. t G J. 



From (HHl) and (lH) we deduce that 



- / iP'u-'d,ig,.^,„ .u)dx+ {^l§£^^ X^^-^)''^'^^^^ ^(T^fyS +^"(^)' 



which in turn implies 

- Jg1p'^U~'^dt{gi-a,n*u)dx V f jrn2, / C*! 



Jgtp'^dx 16r2^jv(S) J SB 



[ {w-W)^dx < ^ +Sn{t), (50) 
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for a. a. t £ J, with some constant Ci = Ci{S, N, v, A) and 5'„(t) = TZn{t)/ /g tp'^dx. 

The fundamental identity (|10l) with = — logy reads (with the spatial variable x being 

suppressed) 

-U~^dt{gi-a,n *u) = -dtigi-a,n * logu) + (log U - l)gi-a,nit) 

+ {- ^Ogu{t - S) + logu(t) + !fc-fhl^) [-9l-a,nM ds. 

In terms oi w = log u this means that 

-U^^dt{gi-a,n *u) = - dt{gi-a,n * w) + {w - l)5l-a,n(0 

+ / 'f{w{t-s)-w{t))[-gi^^^n{s)]ds, (51) 
Jo 

where ^^'(^/) = — 1 — Since 5* is convex, it follows from Jensen's inequality that 

Jg^^'i'(w{t — s, x) — w{t, x)) dx ^ ^ / fg ip^ (w{t — s , x) — w{t, x)) dx 
ip'^dx ~ \ ijj'^dx 

Using this and ([5T|) we obtain 



Jgip'^u '^dt{gi-a,n*u)dx 



Ib "P^dx 

+ / ^{W{t-s)~W{t))[-gi^^^^{s)]ds 
Jo 

^-e-'^dtigi-c^n^e'^), (52) 

where the last equals sign holds again by ([5T|) with u replaced by e^. From ([50)1 and ([5^ we 
conclude that 



/ (w - Wfdx < dt{gi-o.,n *e^)+^ + Suit), a.a. t e J. (53) 

JSB 

lakes sense, s 

C{J;L2{B)) and 



16r2/XAr(B) 
We choose 

. 9/ ^ /■ (54) 
This definition makes sense, since g\-a * ^ G C'(J). The latter is a consequence of g\-a * u S 



- Js^ix)^dx ' a.a.tej, 



where we apply again Jensen's inequality. 
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To prove (HH) and one of the key ideas is to use the inequahties 

PN+i{{{t, x) e : w{t, x) > c{u) + A}) 

< fiN+ii{it,x) G X_ : w{t,x) > c{u) + X and W{t) < c{u) + A/2}) 

+ fiN+ii{{t,x)eK^:W{t)>ciu) + X/2})=:h+l2, A > 0, (55) 
^j,N+i{{{t,x) e K+ : w{t,x) < c{u) - X}) 

< fiN+ii{it,x) e K+ : w{t,x) < c{u) - X and W{t) > c{u) - A/2}) 

+ ^iN+li{it,x)eK+■. Wit) <ciu)~X/2})=:h+h, A > 0, (56) 

and to estimate each of the four terms Ij separately. 

We begin with the estimates for W. To estimate I2 and li we adopt some of the ideas 
developed in ^31'. We set J- := (0,77rr^/"), J+ := {rjTr'^/°' ^Tr^^"), and introduce for A > the 
sets J- (A) ■.= {tC J- : W{t) > c{u) + X} and J+(A) := {t e J+ : W{t) < c{u) - A}. 

Interestingly, positivity and integrability of the function are sufficient to derive the desired 
estimate for /2, cf. also [31] Theorem 2.3]. In fact, with p — rr^/" we have 

e^^ll{J-{X)) = e^pi^{{t e J- : e^(*) > e^'^^e^}) = / dt 

Jj-{X) 

Jj-iX) J J- 



(gi_„*e^)(77p) Jo 

, 32-a(w) 1 r'' , ,x Wit) J. 

[gi-a * )(r]p) 5i-a(w) ./o 



r(l — a) Tyrr 



2/q 



I (2 — a) 1 — a 

h = pi{J-{X/2))fiN{6B) < ' r^'^PN{B), A > 0. (57) 

(1 — a)X 

We come now to 74. For m > define the function Hm on M by Hm{y) — y, y ^ rn, and 
H,n{y) — m + (y — m)/{y — rn+ 1), y > m. Then iJ^ is increasing, concave, and bounded above 
by m + 1. Further, we have Hm G C^(R), and so by concavity 

0<yHU{y)<Hm{y)<m+l, y>0. (58) 

Multiplying ([55)) by e^iJ^(e'^) and employing ([55)) as well as the fundamental identity (ITUl) . 
we infer that 

5t(gi-a,„ *ff™(e^)) + > -^„e^7f:„(e^), a.a. t e J. (59) 

For i e J-|_ we shift the time by setting s = t — rjrr^^"' — t — rjp and put /(s) = f{s + ijp), 
s e (0, (1 — ?y)p), for functions / defined on J+. By the time-shifting identity (IT5)) . (15^1) implies 
that for a.a. s S (0, (1 — ?])p) 

9,(.gi_„,„ * i7,„(e^)) + ^ H™(e^) > T„,„(.) - 5„e^i7;, (e^) , (60) 



and therefore 



21 



with the history term 

T„,m(s) = / [ - gi-a,n(s + ?7p - cr)]i7,„(e^'^°"^) dcr. 
Jo 

For 9 >Q define the kernel rafi G Li^;oc(R+) by means of 

ra,e{t) +0{ra,e * 9a){t) = ga{t), t > 0. 

Observe that r^^o — 9a- Since ga is completely monotone, Va^e enjoys the same property (cf. 
[131 Chap. 5]), in particular ra,e{s) > for all s > 0. Moreover, we have (see e.g. [3T] ) 

r„,e(s) =r(a)5„(s)S„,„(-es"), s > 0, 

where Ea,i3 denotes the generalized Mittag-Lefher-function defined by 

^0 r(na + /3) 

We put 9 = Ci/r^ and convolve (|60l) with Ta^g. We have a.e. in (0, (1 — ry)p) 
ra, e * <9s (^5i^a,„ * iJ„i (e^) ^ = 9s (^q., g * 5i_a,„ * iJm (e^) ^ 

= ds[[ga -6{ra,e * (Ja)] * gi-a,n *ffm(e^)^ 

and so we obtain a.e. in (0, (1 — rj)p) 

+ 9K*r^^e*H^ (e^) - 9r^^ g * (e^) . (61) 
Sending n — oo and selecting an appropriate subsequence, if necessary, it follows that 

-ffm(e'^) >r„,e*T„, a.a. se (0,(l-?7)p), (62) 

where 

f-np 

(s) = / [ - <7i_o (s + 77P - )] i^m (e'^^'O d^. 
Jo 

Observe that for s G (0, (1 — ry)/?) we have 

< 61s" < ^ (1 - ?7)"(Tr2/")" = Ci(l - 77)"t" cj, 

and thus by continuity and strict positivity of Ea^a in (— oo,0], 

ra,e{s)>T{a)ga{s) min ^(-y) =: C2(a, w)r(a)ga(s), s G (0, (1 - ?7)/9). 

S/6[0,a;] 
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We may then argue as in f^P, Section 2.1] to obtain 

gn.(e'^<^0 > C2{a,u) {vpr-'{gi-a*H^{e'^)){f^p), a.a. s G (0, (1 - 77)^)- 

Evidently, Hm{y) y as m -> 00 for all y G M. Thus by sending m — > 00 and applying Fatou's 
lemma we conclude that 

e^(-) > C2{a,u) "I'^^^if^^" M"-i(gi_„ a.a. s G (0, (1 - 77)^). (63) 



We then employ (|63|) to estimate as follows. 

e^Aii(j+(A)) = e^^i({i G J+ : e^(*) < e=(")e-^}) = / dt 

ij+(A) ij+ 



ag2-a{iip) 



ds 

Cola tj)-ifw)^"" /■(i^'')'' 



r(2-a)w " — 



CT "{1 + a) da — C3{a,ri,uj)p. 



Hence 



h ^ Pi{J+{X/2))pn{SB) < 2g3(a,^r,,a;)J ^2/^^^^^^^ ^ > q. (64) 

We come now to h. Set Ji(A) = {t G J- : c - W^(t) + A/2 > 0} and ^^^"(A) = {x G (5S : 
w{t, x) > c + X}, t € Ji{X), where c = c{u) is given by ([M)) . For t G Ji(A), we have 

w(t, a;) - W{t) > c - W^(t) + A > A/2, a; G 

and thus we deduce from (l53l) that a.e. in Ji(A) 

- (c-.y + A)^ (»-"'ab.-.. . «"') + ^ + (65) 

Set x{t,X) = /iAT (rijT (A)) , if t G Ji(A), and x(^i '^) = in case i G J- \ Ji(A). Let further 
i/(y) = {c-logy + Xy^, < y < y* e'^+^Z^. Clearly, H'{y) = {c-logy + Xy^y^^ as well as 

H"{y)=-, — r^-TT^( — r^— T-i)' o<y<2/„ 

(c - logy + A)"'y"' Vc - log y + A / 

which shows that H is concave in (0,y*] whenever A > 4. We will assume this in what follows. 

We next choose a extension H oi H on (0,oo) such that H is concave, < H'{y) < 
H'{y^), y*<y< 2y,, and H'{y) = 0, y > 2y*. Then 

< yS'iy) < p y > 0. (66) 
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In fact, for y £ (O,?/*] we have 

while in case y € 2y^] we may simply estimate 

yH'{y)<2y,H'{y,)< j. 

It is clear that H is bounded above. There holds 

H{y) <j, y>0. (68) 

To see this, note that since H is nondecreasing with H'{y) = for all y > 2y», the claim follows 
if the inequality is valid for all y G [y*, 2y*]. For such y we have by ([S7|) and by concavity of H 

H{y) < H{y,) + H'{y,){y ~ y,) < H{y,) + y,H'{y,) < |. 
Observe also that 

Since H' > 0, and e-^dt{gi^a,n * e^) + C\r-'^ + Sn > on J_ by virtue of dSS]), we infer from 
(|65|) and §^ that 

x(i,A) <e'^i7'(e^)(e-'^at(5i-o.n*e^)+ ^ + 5„) 

< ^'(e^)9,(gi_„,„ * e^) + ^ + ^1^, a.a. t £ J_. (69) 

At a 



16r2/iw(^) 



Since H is concave, the fundamental identity pUj) yields 

i7'(e'^)5t(5i-o,„ * e^) < * ^(e^)) + ( - Hie"^) + ff'(e^)e'^)5i-a,n 

< dt(^gi-a,n * H{e^)^ + jgi-a,n, a.a. t e J-, 
which, together with (|69| . gives a.e. in J_ 

x(M)<a((.._.,...H(e.'')) + II + (70) 

We then integrate ([70|) over J_ = (0, ?7p) and employ ([68]) for the estimate 

* -ff(e^)j(?7p) < - / gi-a,n{t)dt. 
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By sending n oo, this leads to 

Ji(A) 



Hence with C5 = max{4T, C4} we find that 

A<55l!^, A>0. (71) 
A 

It remains to derive the desired estimate for J3. To this purpose we shift again the time by 
putting s — t ~ rjp, and denote the corresponding transformed functions as above by W , w, ... 
and so forth. Set further J+ := (0, (1 — ri)p). By the time-shifting property (ITSI) and by positivity 
of e^, relation ((53)) then implies 

/ I [w- Wfdx < e-^a,(gi_„,„ * e'^) + ^ + Sn{s), a.a. s e J+. (72) 

Next, set J2(A) = {s e J+ : ty(s) -c + A/2 > 0} and f7+(A) = {x e 55 : w(s,x) < c- A}, s £ 
J2(A). For s G J2(A), we have 

W{s)~w{s,x) > W{s)-c + X> A/2, X £ r2+(A), 

and thus ([72]) yields that a.e. in J2(A) 

/ MA^(»+(A)) < - ^ (e-'^g,(gi-.,„.e^)+ ^ (73) 

16r^^Ar(S) (M^-c + A)2V J 

We proceed now similarly as above for the term Ii. Set x(s. A) = /i7v(f^^(A)), if s G J2(A), 
and x(s, A) = in case s G J+ \ Ji(A). We consider this time the convex function H{y) = 
(logy - c + A)~^ for y > y* eC--^/2 ^jth derivative H'{y) — -(logy - c + A)^^y~^ < 0. We 
define a extension H oi H on [0, 00) by means of 



Hiy) 



H'{y*){y - y*) + H{y^) : < y < y* 

H{y) ■y>y*- 



Evidently, ~H is concave in [0, 00) and 

'^-^'^y^y^ (logy.-c + A)2 ^ (W-^' .>0,A>1. (74) 

We will assume A > 1 in the subsequent lines. 
Observe that 

-e^(^)i7'(e^(^)) = , a.a. s G J2(A). 
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Since -H' > 0, and e-^9s(gi-a,« * e^) + Cir-^ + S^^ > on J+ due to d?!]), it thus follows 
from ^ and dTJ) that 



16r2/x^(B) 



< -ff (e'^^j^sigi-a,™ * e"*^) + — + , a.a. s e J+. (75) 

At a 

By concavity of —H, the fundamental identity (fTO|) provides the estimate 
a.e. in J+, which when combined with (|75p leads to 

for a. a. s G J+. We integrate this estimate over Jj^ and send n c» to the result 

J2(A) 



fiN[^J{X))ds = x{s,X)ds< [^jg2-a[(l-V)P) + ) 

-(6.g2-a((l - ?7)Tj +4t;i(l - 77)t) =: Ce r , A > 1. 



Hence with Cy — maxjr, Cg} we obtain that 

,^,Cjf^^ A>0^ (76) 
Finally, combining (l55|), dSS]), and (HZ]), jMl), dUl), dZS]) establishes the theorem. □ 
3.4 The final step 

We are now in position to prove Theorem ll.il Without loss of generality we may assume that 
u > s for some e > 0; otherwise replace u by u + e, which is a supersolution of ([Ij with uq + e 
instead of uq, and eventually let e — > 0+. 

ForO < o- < 1, we set = {tQ + {2~a)Tr^/°' ,tQ+2Tr^^°')xaB and [/; = (<o, io+cr-rr^/") xcrS. 
Clearly, Q_{to,Xo,r) = U'g and Q+(to,xo,r) = Us. 

By Theorem O 



ess sup u 



^ ( _ ^.)ro j 1^ \l,{u^)^ S<a <a<l, J e{0,l]. 



Here C = C(z^, A, 6, r, a, A^) and To — To{a, N). This shows that the first hypothesis of Lemma 
is satisfied by any positive constant multiple of with /3o = oo. 
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Consider now fi=u ^e"^*^"^ where c{u) is the constant from Theorem 13.31 with = U[ and 
= Ui. Since log/i = c(u) — logu, we see from Theorem 13.31 estimate p5|) . that 

fiN+ii{it,x) e Ui : log f lit, x) > A}) < MfiN+iiUi)X-\ A > 0, 

where M = Miiy, A, 6,T,r], a, N). Hence we may apply Lemma 12.31 with /?o = oo to /i and the 
family U^; thereby we obtain 

ess sup /i < All 

Us 

with Ml — Mi^iy, A, (5, r, ?7, a, N). In terms of u this means that 

e^W < Ml essinf u. (77) 

Us 

On the other hand, Theorem 13.21 yields 

Here C — C{i',A,6,T,a,N,p) and ti — Ti{a,N). Thus the first hypothesis of Lemma [2.31 is 
satisfied by any positive constant multiple of u with — p and rj — \/k. Taking /2 — Me""^^"-* 
with c{u) from above, we have log/2 = logu — c{u) and so Theorem 13.31 estimate (jH]), gives 

/iAr+i({(t,x) e U'l : \ogh{t,x) > A}) < MfiN+i{U[)\-\ A > 0, 

where Af is as above. Therefore we may again apply Lemma 12.31 this time to the function /2 
and the sets C/^, and with (3o — p and 77 — we get 

\f2\LAU',)<M2^iN+l{U'lf'^, 

where M2 = M2{v, A, 5, r, 7?, a, Rephrasing then yields 

^iN+l{U'l)-^'^\u\L^^u',) < M2e<-\ (78) 
Finally, we combine ([77|) and ([78l) to the result 

^lN+l{U[y^/P\u\L^^u',) < M1M2 ess^inf M, 
which proves the assertion. □ 

4 Optimality of the exponent 2+Na-2a weak Harnack 

inequality 

In this section we will show that the exponent 2+Na-2a Theorem 11.11 is optimal. 
To this purpose consider the nonhomogeneous fractional diffusion equation on 

d^u - Aw = /, t e (0, T] , X e R^, (79) 
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with initial condition 

u(0,x)=0, a;eR". (80) 

Following [TU], we say that a function u e C{[0,T] x M^) n C((0, T]; C2(R^)) with * m e 
Ci((0, T]; C(R^)) is a classical solution of the problem §^ ii u satisfies ^ and (EOl). For 
any bounded continuous function / that is locally Holder continuous in x, there exists a unique 
classical solution u of the problem (|79p , (l80l) , and it is of the form 

uit,x)= f f Y{t~T,x-y)fiT,y)dydT, (81) 

Jo JR" 

where 

y(t,x)=c(iV)|xr^t"-li/?0(^l^-.|^|2|Ko)^ 

cf. [TU]- Here ^^12 (^l(Ar/2 i) (1 1)) denotes a special H function (also termed Fox's H function), 
see [171 Section 1.12] and [TO] for its definition. It is differentiable for 2: > 0, the asymptotic 
behaviour for z — > cxd and z — > +0, respectively, is described in [101 formulae (3.9) and (3.14)]. 
It has been also proved in [TD| that Y is nonnegative. 

We choose a smooth and nonnegative approximation of unity {(/)„(t, 2;)}neN in R+ x R^ such 
that each (/)„ is bounded. Put / = 0„ in ([75]) and denote the corresponding classical solution of 
([7^ . (ISn)) by Un- Evidently, u„ is nonnegative and satisfies 

d^Un- AUn = (t)n>0, t G (0, T] , X G . 

Hence u„ is a nonnegative supersolution of ([75]) with / = for all n G N. 

Suppose the weak Harnack inequality ^ holds for some p > 2+Na-2a ■ Then, by taking 
Q_ = (0, 1) X B{0, 1) and Q+ = (2, 3) x B{0, 1) it follows that 

(/ uld^lN+lf^''<Cm^Ur,, nGN, (82) 
Jq- <3+ 

where the constant C is independent of n. Since u„ — > y in the distributional sense as n — > cxd, 
we have 

infu„< ^ [ u„d/iAr+i<H ^ ( YdfiN+i <oo, n>no, 

for a sufficiently large tiq. On the other hand, the left-hand side of ([5^ cannot stay bounded, 
since Y (f. Lp{Q.) for p > 2+Na-2a ■ ^^ct, writing Hl^{z) = HI^{z\[nJIi),{i,i)) short, we 
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have 



^ lO "'5(0,1) 

Cl 



>Cil eiN-Np)/2+(^~l)p^^ f p^-^-^PHl^{pyAY dp 



^0 
1 i-t' 



<=■/■/ 

1 



>C2 / r(^-^rt/2+("-l)Pdt, 





with some positive constant C2. The last integral diverges for allp > 2+Na-2a • Hence ([5^ yields 
a contradiction. 

5 Applications of the weak Harnack inequality 

The strong maximum principle for weak subsolutions of ([l} may be easily derived as a conse- 
quence of the weak Harnack inequality. 

Theorem 5.1 Let a G (0, 1), T > 0, and C be a bounded domain. Suppose the as- 
sumptions (H1)-(H3) are satisfied. Let u & Za be a weak subsolution of {Ip in Ot and assume 
that < esssupj^^tt < oo and that esssupQUo < esssupj^^u. Then, if for some cylinder 
Q = {tQ.ta + Tr'^/°') X -6(2:0, r) C VLt with to, r, r > and B{xq, r) (Z fl we have 

ess sup u — ess sup w, (83) 

Q Ot 
the function u is constant on (0,io) x £7. 

Proof: Let M — esssup^^ u. Then w := M — m is a nonnegative weak supersolution of ([T]) with 
uq replaced by uq := M — mq > 0. For any < < + yyr^/" < t^ the weak Harnack inequality 
with p — 1 applied to v yields an estimate of the form 



-{N+2/a) 



/ / {M -u)dxdt<C ess inf (M - u) = 0. 



This shows that u — M a.e. in (0,to) x B{xo,r). As in the classical parabolic case (cf. ^9j) the 
assertion now follows by a chaining argument. □ 

We next apply the weak Harnack inequality to establish continuity at t = for weak solutions. 

Theorem 5.2 Let a G (0, 1), T > 0, and ft C be a bounded domain. Suppose the assump- 
tions (HI) and (H2) are satisfied. Let u E Za be a bounded weak solution of (Qp in fix with 
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uq = 0. Then u is continuous at (0, Xo) for all xq G and ^vccnt x)'^{o.xo) u(t, x) — 0. Moreover, 
letting rj > we have for any cylinder Q{xq, rg) :~ (0, r]rQ^°') x B{xq, tq) C fix cmd r G (0, tq] 

essoscu < c( — ) \u\l^{Qt)^ (84) 

with ess oscq(2.q ,,) = ssssupg^^^^ ^-j — essmfgi-^.^^) and constants C = C(z^, A, 77, a, A^) > and 
S = S{iy,A,Tj,a,N)e{0,l). 

Proof: Let u G Zq, be a bounded weak solution of ([T]) in JIt with uo = 0. Set u(t, x) — and 
x) = Id for t < and a; G i^. For Tb > we shift the time by setting s — t + Tq and put 
f{s) = /(s - To), s G (0,T + To), for functions / defined on (-To,T). Since Du{t, •) = for 
t <0 and 



dt{gi-a.n*u){t,x) ^ dt / 5i_Q.„(t - T)M(r, a;)dr = 9s(gi_a,„ * ?i)(s,a;), 

J -To 

the function tt is a bounded weak solution of 

dfu~dW{A{s,x)Du) =0, s G (0, T + Tq), a; G r?. 

Next, assuming r G (0,ro/2] we introduce the cylinders 

Q,(a;o, r) = ( - r^r^/" .r^r^'^) x 5(^0, r), 

g_(xo, r) = ( - r;(2r)2/", -,y(3r/2)2/") x S(xo, r), 

and denote by Qt:{xo,r) resp. Q-(a;o,'') the corresponding cylinders in the {s,x) coordinate 
system. Let us write Mi — esssupg^^^^ u and mi ~ cssinfg^^^^^ ^^.j u for i = 1,2. Choosing 

^0 ^ 7y(2r)^/", we may apply Theorem 11.11 with p = 1 to the functions M2 — u, u — m2, which 
are nonnegative in (0, 77(2r)^/" + Tq) x B{xo, 2r), thereby obtaining 

/ (M2 - d^w+i < C{Al2 - Ml), 

jQ_(a;o,r) 
/■ (~ _ ^^^^^ < ^(^^ _ 

where C > 1 is a constant independent of u and r. By addition, it follows that 

M2 -m2 < C{M2 - m2 + mi " Ml). 
Writing w(xo,r) = ess supg^^^^^^^-, u ~ essinfg^^^^ .^j u, this yields 

LoixQ, r) < eLo{xQ, 2r), r < ro/2, (85) 
where 6 — 1 — C^^ G (0, 1). Iterating (|85p as in the proof of 11, Lemma 8.23] we obtain 

1 / r N log9/log(l/2) 

w(a;o,»-) < T I — '^(a;o,»-o), J' < ^-o- 
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The estimate (IMl) then fohows by transforming back to the function u and using that m = for 
negative times. In particular, we also see that u is continuous at (0,xo) for all xq € n and that 

lim(^t,x)^(o,xo) u{t,x) = 0. □ 

The last application is a theorem of Liouville type. We say that a function u on IR+ X is 
a global weak solution of 

d^u - div {A{t, x)Du) = 0, (86) 
if it is a weak solution of ^ in (0, T) x B(0, r) for aU T > and r > 0. 

Corollary 5.1 Let a e (0,1). Assume that A G ico(R+ x M^;M^^^) and that there exists 
ly > such that 

{A{t,x)(,\^) > i^l^l"^, for a.a. {t,x) 6 R+ x M^, and all (, € M^. 
Suppose that u is a global bounded weak solution of i86]) . Then u — a.e. on x M^. 
Proof: For r > and xo = it follows from the proof of Theorem 15.21 that 

w(0,r) < 6la;(0,2r), r > 0, (87) 
where 6 £ (0, 1) is independent of r and u. By induction, (|87|) yields 

^(0,r) < rw(0,2"r) < 2^ r > 0, n e N. 
Sending n oo shows that u is constant. The claim then follows by Theorem 15.21 □ 
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